Oscillatory Instability in Two-Dimensional Dynamic Fracture 
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The stability of a rapid dynamic crack in a two dimensional infinite strip is studied in the frame- 
work of Linear Elasticity Fracture Mechanics supplemented with a modified principle of local sym- 
metry. It is predicted that a single crack becomes unstable by a finite wavelength oscillatory mode 
at a velocity Vc, O.Scr < Vc < 0.85c_r, where cr is the Rayleigh wave speed. The relevance of this 
theoretical calculation to the oscillatory instability reported in the companion experimental Letter 
is discussed. 



Introduction: High precision experiments on dy- 
namic fracture in slabs of amorphous materials re- 
vealed very interesting instabilities in the form of micro- 
branching. Above a critical velocity of about OAcr, 
where c_r is the Rayleigh wave speed, a single, straight, 
rapidly moving crack is unstable against the appearance 
of small diameter side-branches that affect both the mor- 
phology and the velocity of propagation 0. An im- 
portant observation regarding this instability is that al- 
though experiments are typically performed on quasi- 
two-dimensional samples (i.e. samples for which the 
third dimension is significantly smaller than the other 
two dimensions), the instability is intrinsically three- 
dimensional 0) Q ; at the onset of instability the micro- 
branches occupy only a small fraction of the small third 
dimension, barely misting the mirror quality of the main 
crack. 

Nevertheless, theoretically there were a number of at- 
tempts to explain this instability in the context of Lin- 
ear Elasticity Fracture Mechanics (LEFM) in two dimen- 
sions, where the experimentally observed instability was 
interpreted as a macroscopic crack bifurcation. It was 
shown that above some critical velocity such a bifurca- 
tion is allowed on energetic grounds 4], but not neces- 
sarily realized dynamically. The aim of this Letter is to 
show that two-dimensional LEFM predicts the existence 
of a dynamical oscillatory instability. We refer the reader 
to the experimental companion paper where a similar 
conclusion is offered on the basis of crack propagation in 
thin films loaded with a fixed grip . Here we construct 
a theoretical analog: a semi-infinite straight crack prop- 
agating at a constant velocity in an infinitely long two- 
dimensional strip under fixed grip boundary conditions. 
The standard framework of LEFM is supplemented with 
a modified principle of local symmetry 6, 7], and see 
below for details. The analysis is based on a recently de- 
rived solution for the linear perturbation problem of the 
dynamic stress intensity factors using the weight func- 
tions method We find that an oscillatory mode 
of finite wavelength becomes unstable above a critical 
velocity Vc, O.Scfl < Vc < O.SbcR. 

The perturbation problem: Consider a semi- 
infinite straight crack dynamically propagating at a con- 
stant velocity v in an infinitely long two-dimensional strip 



of width 2h under fixed-grip boundary conditions. A co- 
ordinate system {x, y) is located on the central line of 
the strip with x being the direction of crack propagation. 
The unperturbed crack configuration Mq at any time t is 
described by 



Afo = {{xq, yo) : -oo < xq < vt, j/q = 0} 



(1) 



We then consider a configuration Af^ that results from a 
small time-dependent perturbation of the crack path 

Me = {{xe,ye) ■ -OO < x^ < vt + eip{t), y^ eip{x)} , 

where ip(t) and tpix) are smooth dimensionless functions 
that define the longitudinal and transverse perturbations, 
respectively. The dimensional amplitude e satisfies < 
e/h 1, ensuring that both the speed and the path are 
only slightly perturbed. For convenience we define a co- 
moving frame by the transformation X = x — vt and 
Y = y. The material is assumed to be linear elastic and 
isotropic, with the displacement field u(a:, y) satisfying 
Navier's equation 



(A + ^J.)V (V • u) + ^V^u = pii 



(3) 



where p is the density and the dots stand for partial time 
derivatives. A and p, are the Lame coefficients that re- 
late the spatial derivatives of u to the components of the 
stress tensor cr 11] 



A(5y ^ dkUk + /i (diUj + djUi) 



(4) 



with i,j = x,y. The boundary conditions are given by 

(JijUj = on the crack, u{x, ±h) = ±wy , (5) 

where n is an outward unit vector normal to the crack 
faces, y is a unit vector in the y direction and w is a 
constant. As the fracture process is localized near the 
crack tip we are mainly interested in the fields in that 
region. To that aim we define a shifted and rotated frame 
{X', Y') according to 



X' + iY' = {X - etp{t) + i{Y- e'tP{vt))}e'^ 



(6) 



where 9 is the angle between the crack tip orientation 
and the X-axis such that tan(6') = edxip{e(p{t) + Vt). 
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Then, the asymptotic expansion of the stress field along 
the Y' = direction, sufficiently close to the crack tip 
{X' 0+) yields [13 

a^,^, {X , 0, t) ~ J^^ ^ ^x'Y' : 0. ^) 



+ 



^{[uj^{Pi'^)-{u,)^[Pi%\x=oAn) 



(7) 

where and K-^^ are the mode I (tensile) and II (shear) 
"stress intensity factors" respectively. Our goal is to ob- 
tain ifj and /iTjj as functionals of the perturbation func- 
tions ip{t) and ip{x) to linear order in e. 

Solution: As a result of the long range nature of elas- 
tic interactions we expect K^^^ and K^p to depend on 
the whole zeroth order solution. Therefore, we write the 
displacement field as u = u^^'>+u^^'> and the stress field as 
(T = CT^'^^ +(T^^\ Then, we represent the unperturbed dis- 
placement field as =u'''^^+ywy /h. Here u^^") is a field 
satisfying Eq. in a strip with a straight crack and the 
boundary conditions (Jxy — O, (7yy — — {2fj, + X)w/h on the 
crack and it'*'' (x, ±/i) — 0. In addition, the asymptotic 
expansion of the unperturbed stress field is given by 



.(0) 



iX,0) 



Af^iv)Vx 



(8) 



where we have included the sub-leading contribution pro- 
portional to ^["^ . Accordingly, the expansion for the dy- 
namic stress intensity factors is written as 

(9) 

Note that there is no zeroth order contribution to K-^-^ 
since the unperturbed crack is pure mode I. 

The zeroth order problem is rather straightforward and 
the resulting solution is presented in The goal now 
is to obtain K^^'' and K^p as linear functionals of the 
perturbation functions tp{t) and ipi^)- The crucial step 
is the construction of an auxiliary field U, whose compo- 
nents are the so-called dynamic weight functions U 
satisfies Eq. in a strip with a moving straight crack 
along the X-axis and possess the following properties: (i) 
[U] (X < 0,0) = and U{X,±h) = 0. (ii) The com- 
ponents of the stress tensor a^^ (U) are continuous and 
a„ (U) {X > 0,0) = 0. (iii) \u^] (X,0) ~ c^X-^/^6{t) 
for X ^ 0^. Here q are constants and [•] denotes the 
jump of a function across Y — 0. The field U is found by 
solving the associated Wiener-Hopf type equations and 
the explicit Fourier space representation of [U] is given 
in lid. With these objects at hand the original reasoning 
of ja ISl applies without change and the linear functional 
ifti) and are shown to be (see Ref. 10] for details) 

K('\v,t) = 6if("){Q,(t)*^(t)} + ey|A(°)^(t) (10) 



where (•) denotes the average of a function across F = 
and * denotes a convolution with respect to all relevant 
variables. The function {U^) can be expressed in terms 
of the components of \U] , while the Fourier transform of 
the field Q is obtained by a two-term asymptotic repre- 
sentation of the Fourier transform of [[/], see 10]. The 
functions Q{v) and T(w) are known universal function of 
V 1^; [P^i)] and {P^^) are given by 



r^{x)pv^[dl,u^:\X))H{-X) , 
{^(x)T(i)(X) - 9, (^(x)T(0)(X))}i/(-X), 



[Pi'' 
T^°\X 



Xw 



<^^(X,0) + — , 



(12) 



Here H{-) is the Heavisidc function and recall that x = 
X + vt. 

The various terms on the RHS of Eq. (|11|) have clear 
physical meanings. The first term is a result of shifting 
the crack tip out of the symmetry line y = to y = iplvt) 
in the presence of the unperturbed asymptotic field of 
Eq. By dimensional analysis it is obvious that this 
term must be proportional to the sub-leading contribu- 
tion A[°'> . The second term on the RHS of Eq. ((TTl) is a 
result of changing the crack tip orientation by an angle 
9 ~ edxi^ivt) (see Eq. ^) in the presence of the un- 
perturbed asymptotic field of Eq. JSj) . These two terms 
are both local and instantaneous. As a result of long 
range elastic interactions and wave propagation we ex- 
pect also non-local contributions in time and space. The 
third term on the RHS of Eq. (|ll|l is a result of the in- 
teraction of the tip with waves that were transmitted by 
itself at earlier times. The precise nature of the interac- 
tion is carried by the time-dependent kernel (^) ■ The 
last term on the RHS of Eq. Ijllll represents the inter- 
action of the tip with all the history of propagation of 
the crack. This interaction is non-local both in time and 
space. The functions [P^^^] and (P^^)) are the efi^ective 
mode II (shear) loads on a straight crack, loads that are 
induced by the broken up-down symmetry due to ip{x). 
With this interpretation in mind, one observes that the 
dynamic weight functions U are response functions that 
quantify the effect of forces applied to the crack line on 
the asymptotic stress field near the tip. An instability, 
if exists, should emerge as result of the competition be- 
tween these various terms. 

Stability analysis: In order to study the stability of 
the straight crack against small perturbations, one must 
have an equation of motion for the crack tip. Such an 
equation does not emerge from LEFM. This theory can 
determine the stress and displacement fields, under the 
linear stress-strain assumption of Eq. Q , given the crack 
and the loading histories, but cannot tell where and at 
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FIG. 1: (color online) hK^p /eK^'^^ as a function of rjh for two 
different velocities. The squares correspond to v = O.Scr and 
the circles correspond to u = 0.85c_r. The onset of instability 
occurs at v^, O.Sch < Vc < 0.85cij, where K^^^ first changes 
sign from positive to negative. The wave number at the onset 
of instability is approximately given by r]ch^lA. 



what rate the crack tip will propagate. The rate of crack 
propagation is determined by assuming that the energy 
flowing to the crack tip per unit crack extension per unit 
time, J, equals vT{v), where T{v) is the "fracture en- 
ergy" . The quantity J is calculated using the dynamic 
stress intensify factors in the asymptotic expansion Eq. 
l|7|l [l^ . Considering a pure longitudinal perturbation, 
■0(2;) = and 1^9(^)7^0, one obtains, to first order in e, 



2T{v) 2V{v) 



(13) 



where Eq. was used and J-{v) is a known universal 
function 10]. In 10] the case of constant fracture en- 
ergy, d^Tiv) = 0, was studied. Since for the fixed-grip 
strip configuration J = vG, where the "energy release 
rate" G is constant, every velocity v is possible as long 
as G = r. It was shown in ^3] that the crack is un- 
stable for all the velocities below v « O.Gcfl. On the 
other hand, real materials have dyT{v) > such that a 
unique velocity is selected according to G — T{v) and 
longitudinal stability is expected 13]. Wc verified that 
Eq. (|13|l indeed predicts longitudinal stability for all ve- 
locities as long as 9„r(u) > 0. Accordingly, we can safely 
assume a straight crack propagation at a constant veloc- 
ity and focus on transverse perturbations, i.e. ipi^) 7^0 
and (p{t) — 0. Note that all the calculations, here and 
below, were performed for a material with a fixed ra- 
tio between A and ^ such that Poisson's ratio for plane 
stress condition was ly = 0.23 |0. The results depend 
smoothly (and weakly) on 1^. 

To study the transverse stability of the straight crack 



80 
di 



(14) 



where f{K-^,K'^^,v) is a positive definite dynamic mate- 
rial function that quantifies how the asymmetry in the 
asymptotic fields Eq. Q , characterized by , is trans- 
ferred to the crack tip itself; is defined in Eq. Eq. 
I|14|l was applied quite successfully in [3,^3]- The equa- 
tion was obtained by using symmetry arguments and as- 
suming (i) that the crack tip is not appreciably blunted 
and the crack path is a trace of a point, (ii) that the 
region around the tip in which linear elasticity theory 
is not applicable (the so-called "process zone") is small 
compared to any other length scale, (iii) the crack path 
is smooth. By dimensional considerations we define a 
dynamic length scale as 



l{v) 



f{K„Kl,v)K, 



(15) 



As this length scale is expected to be small, in our case 
implying i[v)lh <C 1, Eq. (|14|) is approximated by many 
authors by K^^jK-^ « 0, which is the well-known "prin- 
ciple of local symmetry" [13. Eq. IHI) is thus a mod- 
ified principle of local symmetry. Note that for £(u) to 
be well defined in the limit of small w, we must have 
/(ifj, i^j^j, f ) C!C V in that limit. 

To proceed we can consider two types of transverse 
perturbations, one localized right at the tip |lOl ]. and the 
other global, like a small amplitude sine function [ifij ]. 
We follow the latter approach since real cracks are never 
perfectly straight. Analyzing the perturbation in linear 
modes, we consider 



i}{X + vt) = sin{r/(A: + vt)} . 



(16) 



Then, we note that Eq. l(T^ . with dt9 — evdxx'>P{vt), is a 
linear integro-differential equation for ip{x). In principle, 
such an equation can be analyzed by standard methods, 
but we encounter a difficulty as f{K^,K-^^,v) is not yet 
known. To overcome this difficulty, we focus, without loss 
of generality, on times such that tplvt) ~ 0. The dynam- 
ical equation (|14|l then predicts that if 9 and K^^ have 
the same sign, with / > 0, then dtO has the opposite sign 
and \9\ decreases, which means that a small perturbation 
decays. By the same argument, for 9 and K^^ having the 
opposite sign a small perturbation grows. This criterion 
is identical to the one used in Refs. Therefore, 
considering the perturbation introduced in Eq. H16(l . we 
interpret a change of K^^ from positive to negative as 
an instability. Note that this criterion is independent 
of the exact form of the function f{K-^,K^^,v) and the 
wavelength of an unstable mode is expected to be cor- 
rect up to modifications of the order of £{vc)- Having an 
explicit instability criterion at hand, we calculate numer- 
ically i^fi) in Eq. ^ as a function of the wave number 
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r] for various velocities and look for the lowest velocity 
for which an intersection with the x-axis occurs. In Fig. 
^we show hK^p / eKj-^^ as a function of rjh for two rep- 
resentative velocities, one below the instability and one 
above. We observe that for v = O.Scr the curve is al- 
ways positive, implying stability, while for z; = 0.85cfl the 
curve changes sign from positive to negative, implying 
instability. The normalized wave number rjch at thresh- 
old is approximately rjch « 1.4. The instability should 
be understood as resulting from a competition between 
local stabilizing terms and non-local destabilizing terms 
in Eq. ifTT)) (see discussion above), which we believe to 
be the dynamic counterpart of the phenomenon studied 
in 

Discussion and summary: While oscillatory insta- 
bilities were observed and discussed before in a number 
of fracture context s in the context of quasi-static 

thermal cracking, [l8j in the context of lattice simula- 
tions, [iflj l in large biaxial strain experiment in rubber, 
[ifij l in context of a phase field model), it appears that 
the present is the first calculation predicting such a high 
velocity instability, purely on the basis of LEFM. This 
prediction is in correspondence with the experimental ob- 
servations in quasi-two-dimensional thin layers where the 
side branching instability is suppressed (see the compan- 
ion Letter (^]). 

Notwithstanding the correspondence of the predicted 
onset of oscillatory instability compared to the exper- 
iment, one should point out that the predicted wave- 
length of the saturated post-instability solution differs 
qualitatively in theory and experiment. The theory, be- 
ing based on LEFM that contains no intrinsic scales, is 
bound to predict wavelength that scales with the width 
h of the strip. The experiment predicts saturated oscil- 
latory cracks with wavelength that do not depend on the 
width of the strip. Interestingly enough, the experiment 
also reports an excellent agreement with LEFM for all 
the properties of the dynamic crack up to the instability 
point. It is tempting to interpret these findings by stating 
that LEFM supplemented with the modified principle of 
local symmetry is trustworthy up to the instability, but 
then fails once the oscillatory motion is stabilized. The 
indication is that while the crack is straight all the non- 
elastic aspects can be lumped into an effective fracture 
energy T(v). The true free boundary dynamics of the 
crack, once deviating from the straight line, call for a 
new understanding that must go beyond LEFM with or 
without principles of local symmetry. It is very possible 
that the actual shape of the tip, including its radius of 
curvature, the distortion with respect to the symmetry 
axes and all other dynamic degrees of freedom become es- 
sential in determining the actual path of the crack. Such 
degree of freedom can introduce intrinsic length and/or 



time scales. The only way to reach such a theory is to 
derive equations directly for the free boundary that de- 
fines the crack. A theory that advances in this direction 
will be presented in due course |2ll |. 
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